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1. Introduction 

The dusty envelope of the carbon star IRC+10216 is a rich radio source where at least 
70 molecules have been observed (see, e.g., Cernicharo et al. 2000). Particularly rich in long 
linear carbon chains, it is the source where the polyacetylenic radicals C„H {n = 4,6,8; see 
Guelin et al., 1987; Cernicharo & Guelin 1996) and C„N (n=3,5; see Guelin et al., 1998) 
were first detected and identified. The butadiynyl radical C4H has been observed there in 
rotational lines of the ground and several vibrationally-excited states. The excited state lines 
are strong, owing to perturbation by the low-lying ^4^11 electronic state (Yamamoto et al. 
1987). Notably poor in cations, IRC-l-10216 has turned out to be the richest source of anions 
in the sky. All known interstellar anions, CgH^, C4H^, CgH", and CsN" (Thaddeus et al. 
2008 and references therein), have been identified there, where they are fairly abundant. 
The presence of carbon chain negative ions in space was predicted long ago (see e.g. Sarre 
1980, Herbst 1981) on the ground that electron radiative attachment is efficient for molecules 
with large electron affinities and dense vibrational spectra. The presence of C3N and C5N 
in IRC+10216 and the recent detection of CaN" in this source, imply that CsN" may be 
present. 



2. Observations 

The astronomical observations presented here are mainly from a spectral line survey of 
IRC+10216 done with the IRAM 30-m telescope between 1995 and 2008 (Cernicharo et al. 
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2000; 2008 in preparation). The 3-mm part of this survey covers the 80-115 GHz band with 
an r.m.s. noise of 0.3-2 mK per 1 MHz-wide channeL Data were taken with the secondary 
nutating by 90" — 120", an arrangement which produced very flat spectral baselines. Two 
SIS 3-mm receivers, with orthogonal polarizations and system temperatures of 100 — 130 K, 
were used simultaneously. Pointing and focus were regularly checked on planets and on the 
strong nearby quasar OJ 287. In Fig. 1 & 2 the intensities are given in Ta*-, the antenna 
temperature corrected for atmospheric absorption and rear spillover losses, by means of the 
ATM code (Cernicharo 1985). A few complementary observations at 22 GHz were carried 
out in March 2008 with the MPIfR Effelsberg telescope, both on IRC+10216 and on TMC-1. 
Derived line parameters are given in Table 1. 



3. Results: Vibrationally excited CgH 

During the survey, about 1500 lines in the 80-115 GHz band were observed, most from 
known molecules. A few hundred escaped identification; among these, the four harmonic 
series with either integer or half-integer quantum numbers in Fig [1] and Fig [2] stand out; for 
one series, B1394, the lines are split into doublets (see Tables 1 & 2). Since the B and D 
values are very close to those of CgH and C5N (B= 1391.2 MHz, D= 51 Hz, and B= 1403.1 
MHz, D= 50 Hz, respectively), the carriers of the four series must be linear molecules of 
weight and size similar to that of these two radicals. 

CgH is a linear radical with a ^H ground electronic state (Cernicharo et al., 1987). It 
has 6 stretching modes, ui — vq and 5 bending modes, z/y — Vu. The modes viq and vu have 
the lowest energies, ~210 and ~93 cm~^, respectively, according to ab-initio calculations by 
Cao and Peyerimhoff (2001); the next mode is at about 400 cm~^. 

Even though excited vibrational states of long carbon chains are susceptible to infrared 
fluorescent excitation in IRC+10216, vibrational temperatures remain low: e.g., 40 — 50 K 
for the z/7 mode of HC3N and the ph mode of HC5N (Cernicharo et al. 1987, 2008). Only 
the lowest bending mode of CgH {vn) is therefore expected to be sufficiently populated to 
detect. The first excited electronic state of CgH, ^S, is known to lie very close in energy to 
the ^H ground state (Murakami et al., 1987). Owing to coupling between the vibrational 
bending modes and the electronic orbital motion of the molecule (the Renner- Teller effect), 
each bending mode of CgH is split into one (with half-integer quantum numbers) and two 
vibronic states (^Ti^ and with integer quantum numbers). According to Herzberg 
(1966), the splitting between the S"*" and S~ states, A^, produced by the Renner- Teller effect 
is 

(1) K^={P^lo+^eWYl\ 
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Fig. 1. — Spectra of IRC +10216 observed with the IRAM 30-m telescope, showing hnes 
from the B1389 series assigned here to CsN". The marginal weak line U83278 is worth 
noting, because it is within 0.1 MHz of the J=l-0 line of CCH^ (see text). 
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Fig. 2. — Spectra of IRC +10216 observed with the IRAM 30-m telescope showing selected 
hnes pertaining to the series B1390, B1394 and B1401. These three series of hnes are assigned 
to vibronic states of the uu bending mode of CeH. 
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where Aso is the spin-orbit constant of the ground ^11 electronic state of CqH (A5'o=-15.04 
cm~^, Linnartz et al., 1999), e is the Renner parameter, and uj is the frequency of the bending 
mode. The and ^S"^ states he above and below the state by A2/2. The values of B 
and 7 for the states, the effective spin-orbit constant for the vibronic state, and the 
effective rotational constants for the two ladders of this state are given by (Herzberg, 1966), 

(2) Bi=B{l ± A|oB/A|), 

(3) 7s=2B(l-2ea;/As ± A|oB/A|), 

(4) Ae//=A5o(l-3eV4), 

(5) Be//A=B(1±B/A,^;), 

where B is the unperturbed rotational constant of the bending mode. These expressions 
have to be taken with caution as additional ro-vibronic coupling with higher lying vibronic 
states may introduce significant higher-order corrections (see, e.g.. Brown, 2003). 

In view of the close agreement between the rotational constant of CeH and those of our 
four new series, it was probable that one or more may by produced by the low-lying i/n bend 
of CeH. By analogy with HC5N, which has a vibrational spectrum similar to that here, CeH 
in its uii bending mode is expected to have a rotational constant 1.004 times larger than 
that in its vibrational ground state, i.e. close to 1395 MHz. The assignment of B1390 and 
B1401, both with half-integer quantum numbers and whose average B value is close to 1395 
MHz, to the ^A3/2 and ladders of the ^A vibronic state of z/n was therefore logical. 

As discussed below, these assignments have been confirmed in the laboratory (see below). 

Because B1394 shows well resolved doublets with integer rotational quantum numbers 
(see Fig [2]), this series is almost certainly from a vibronic state. The most likely carrier 
is the lowest of the two ladders of the i^u state, the ladder. We note that we have 
found no other series of line doublets of comparable or greater intensity that could match 
this ladder, would B1394 pertain to another vibrational state. 

With the assignment of B1394 and B1390/B1401 series to the ^S" and the ^A ladders of 
uii, we can derive from relations (1) through (5) a set of ro- vibrational parameters that give 
a good fit to the data: Ba=1395.72 MHz, Ae//^=11.13 cm-\ £^=0.352, As ~145.0 cm'^ 
and uj ~120 cm~^, a value compatible with that calculated for uu by Cao and Peyerimhoff 
(2001). Therefore, under the assumption that relations (1) through (5) apply to CgH, the 
^A, and ^S"*" vibronic states lie 68, 173 and 277 K above the ground state. The ^S"*" 
vibronic state will have a rotation constant (B+) so close to that of the state (B~), and 
a doublet separation so similar to (7+ ~ 7") that its rotational lines wiU be blended 
in IRC-l-10216 with those of up to very high J (where lines are probably too weak to 
detect). The column densities derived for the z/n vibronic states are: N(^S~)~1.4xN(^A)=8 
10^^ cm~^. Both vibronic states have a rotational temperature of ~45 K — somewhat higher 
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than that measured for the ground vibrational state (31±2 K; Cernicharo et al., 2007). The 
column density in the ground state is calculated to be 6.6 x 10^^ cm~^, indicating that about 
20% of CeH is in the uu state. 

Recent laboratory measurements at both millimeter- and centimeter-wavelengths con- 
firm that B1390, B1401, and B1394 arise from and components of a low-lying vibra- 
tionally excited state of C^H (Gottlieb et al. 2006). The component has been observed 
over a wide range of rotational excitation (A^=5 to 85), under experimental conditions which 
optimize lines of ground state CqH. To reproduce the observed spectrum, an effective Hamil- 
tonian with several higher-order terms in the spin-rotation constant is required. Although 
weaker by a factor of three or more relative to B1394, the high-lying component has 
also been observed in the same laboratory discharge. The spectroscopic constants derived 
from the astronomical data agree to within la with those derived from the larger and more 
precise set of laboratory data for both the ^A and ^E components. A complete account of 
the laboratory observations will be given elsewhere. 

There is so far no laboratory evidence for the remaining astronomical series B1389 under 
conditions where lines of ground state and vibrationally excited CgH or ground state CgH" 
are observed. 

4. Assignment of B1389 with CgN" 

The remaining series, B1389, has integer J and exhibits no evidence of hyperfine or other 
structure (see the J=35 — 34 hue in Fig. 1). For this reason, and because of the strength of 
these hues, it cannot be assigned to any of the bending modes of CeH. Its carrier is almost 
certainly a new linear molecule with a ^E ground electronic state, and a molecular weight 
close to that of CeH and C5N. 

The most obvious candidates are the ions CgH", C6H+ and CsN". CeH" has already 
been detected in the laboratory and in IRC-l-10216 (McCarthy et al., 2006), and is immedi- 
ately ruled out. Linear C6H+ has a ^E ground electronic state (according to Fehrer & Maier, 
1994), and only a fraction of its rotational transitions will be harmonically related with 
B~1340 MHz. Moreover, only one cation, HCO+, has so far been detected in lRC+10216, 
and it has a very low abundance. We therefore conclude that C6H+ is not the carrier of 
B1389. 

Finally, CsN^ has the right ground electronic state and calculated rotational and 
distortion constants: 1389 MHz and 33 Hz (Botschwina, 2008, private communication; see 
also Aoki 2000). Botschwina quotes an error as small as .5% for the rotational constant. 
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The match with B1389 is perfect, and it is therefore tempting to conclude that C5N is our 
new molecule. 

The only reservation we might see to this identification is the intensity of the lines, which 
in IRC+10216 are about twice stronger than those of the parent species C5N. Supporting 
the identification, the lines of the anion benefit from a more favorable partition function (no 
doublets) and from a larger permanent dipole moment: the latter has been calculated by 
Botschwina (2008) to be 5.2 D. Assuming our identification is correct, we derive a rotation 
temperature of 37±6 K and a C5N~ column density of 3.4x10^^ cm~^ (see Fig [3]). Neutral 
C5N has been observed in IRC+10216 by Guelin et al. (1998), who derived a column density 
of 6x10^^ cm~^ for a dipole moment of 3.39 D (Botschwina 1996). The abundance ratio 
between the neutral and the anion is then only 1.8, the largest relative abundance observed 
so far for any anion. 

A number of clues, however, suggest that we may have underestimated the abundance 
of neutral C5N. First of all, as noted by Guelin et al. (1998), the abundance decrement 
C3N/C5N is an order of magnitude larger than that of HC3N/HC5N and C4H/C6H, both in 
IRC+10216 and in TMC-1 — yet all chemical models predict it to be similar. Second, in 
the laboratory the lines of C5N are far weaker than those of CeH, yet the production of both 
should be similar. Pierre Valiron (private communication 2008) has re-calculated the ground 
state symmetry and dipole moment of C5N at restricted open-shell Hartree-Fock and coupled 
cluster level, using augmented correlation consistent Dunning's basis sets. Using a double 
zeta set, he derives a ^11 ground state with a small (~ 1 D) dipole moment, in agreement with 
the previous calculations of Pauzat et al. (1991); alternatively with a quadruple zeta set, he 
finds a ground state and a large dipole moment (~ 3.4 D), in agreement with Botschwina 
(1996). Valiron concludes that the dipole moment of C5N may well lie between those two 
values in the case of admixing between the and ^11 states similar to those observed for 
C4H and CgH. A dipole moment twice smaller than that calculated for the unperturbed 
state would raise the CsN/CsN" abundance ratio to 8, making it only slightly smaller than 
the CeH/CeH" ratio (see Cernicharo et al. 2007). 

To further investigate the formation of C5N and CsN", we have modeled the chemistry in 
the external layers of the circumstellar envelope of IRC+10216 with the same time-dependent 
chemical model we have used to calculate the abundances of C4H and C4H~ (Cernicharo et 
al., 2007) and C3N and C3N~(Thaddeus et al., 2008). Fig H] shows the computed abundance 
as functions of the distance to the star, R. In these models we have assumed an electron 
radiative attachment rate for C5N of 2xl0~'^x (T/300)~*''^ cm^ s~^, close to those adopted 
for the largest CnH radicals. In the region of interest (i? = 3 — 5 x 10^^ cm), C5N is found to 
be 2-10 times more abundant than CsN", which is in good agreement with our identification 
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of B1389 with CsN". We note that our spectral survey should have been sensitive enough 
to detect another series with a rotational constant close to that predicted by Botschwina, 
assuming that B1389 does not arise from this anion and that the dipole moment of C5N~ is 
as large as 5.2 D. No such series was found. 

Laboratory searches to detect CsN" have also been undertaken with the microwave 
spectrometer at Harvard on the assumption that this anion is the carrier of B1389. No 
frequency search is required, because the rotational constant derived from the astronomi- 
cal data is good to a few kHz, and the nitrogen hfs is negligible throughout much of the 
ccntimcter-wave band. As for CeH", experimental parameters were first optimized to pro- 
duce strong lines of the corresponding neutral radical, C5N. A number of searches for CsN" 
were then performed using different discharge voltages (ranging from 1000 V to 600 V); still 
other searches using a range of concentrations and gas pulses of various lengths, were also 
attempted, all without success. Our failure to detect C5N~ is not particularly disturbing, 
because finding C4H~ and CaN" was very difficult (McCarthy & Thaddeus 2008), requir- 
ing repeated trial and error. Additional experiments are needed to find the highly specific 
conditions which are apparently required to make this anion. 

Finally, a number of other candidates for B1389 were considered and rejected, including 
AlCCCN, AINCCC, NaCCCN, and NaNCCC. All have electronic states, large dipole 
moments (/ia=6.2, 4.0, 13.6, & 16.6 D, repectively), and about the right mass and size. The 
calculated rotational constants, however, differ by more than 4% from 1389 MHz, versus an 
estimated uncertainty of 1% (Pctric 1999) — probably a prohibitive discrepancy. Finally, 
MgCCCN, whose calculated rotational constant is close to B1389 (1397 MHz - see Petrie 
1999), has a not a ground electronic state. 

With no plausible alternative, we therefore conclude that B1389 is almost certainly 
CsN". The observation of low- J lines in the cold source TMC-1 would strengthen this con- 
clusion. We have searched for the J=8-7 transition in TMC-1 with the Effelsberg telescope, 
and have detected a 3a feature at the right frequency; a further factor of two improvement 
in sensitivity is required for confirmation. The final proof, however, must await detection of 
the B1389 lines in the laboratory. 

Among the unidentified lines of Fig. 1, we note that U83278 is within 0.1 MHz of the 
J=l-0 line of CCH~. Assuming the line arises from this anion, we calculate an abundance 
ratio CCH/CCH- ~12500. The J=3-2 CCH" transition would then be strong enough to be 
detectable; a search for this transition to confirm the present assignment is under way. 

Shortly before this Letter was submitted, Pierre Valiron kindly calculated for us the 
structure of C5N", and gave an estimate of its dipole moment. This was one of his last 



- 9 - 




5 10 radius (") 30 50 




10 radius (cui) 10 

Fig. 4. — Model abundances of the neutral radicals C^H, C„N and their anions in the outer 
envelope of IRC+10216. 
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scientific contributions; he died untimely on August 31*'*. Pierre was a valued colleague 
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Table 1. Observed line parameters in IRC +10216 







Obs. Preq. 


Cal. Preq. 




Vexp 


Transition 


(MHz) 


(MHz) 


(K km/s) 


(km/s) 


B1389 (C5N-) 


J= 


=29-28 


80550.7(4) 


80550.7 


0.30(3) 


14.6(8) 


J= 


=30-29 


83328.0(3) 


83328.0 


0.23(3) 


15.3(7) 


J= 


=32-31 


88883.1(10) 


88882.7 


0.20(3) 


14.5" 


J= 


=33-32 


91660.0(6) 


91660.0 


0.25(3) 


15.4(7) 


J= 


=34^33 


94437.6(5) 


94437.3 


0.35(3) 


15.6(8) 


J= 


=35-34 


97214.7(3) 


97214.5 


0.21(5) 


14.8(6) 


J= 


=36-35 


99992.3(8) 


99991.8 


0.22(3) 


15.1(8) 


J= 


=37-36 


102769.5(10) 


102769.0 


0.24(4) 


14.5" 


J= 


=38-37 


105546.4(10) 


105546.1 


0.17(3) 


14.5" 


J= 


=39-38 


108323.0(5) 


108323.2 


0.15(2) 


14.5" 


J= 


=40-39 


111099.1(15) 


111100.4 


0.16(3) 


14.5" 


B1390 


J= 


=59/2-57/2 


82007.0(5) 


82006.4 


0.29(2) 


16.8(9) 


J= 


=61/2-59/2 


84786.4(2) 


84786.5 


0.23(2) 


13.7(7) 


J= 


=65/2-63/2 


90346.9(10) 


90346.9 


0.14(2) 


14.5" 


J= 


=67/2-65/2 


93127.0(10) 


93127.1 


0.15(5) 


14.5" 


J= 


=69/2-67/2 


95907.3(4) 


95907.3 


0.24(2) 


14.5" 


J= 


=73/2-71/2 


101468.5(10) 


101467.9 


0.23(4) 


14.5" 


J= 


=75/2-73/2 


104248.2(10) 


104248.2 


0.17(2) 


14.5" 


J= 


=77/2-75/2 


107028.7(5) 


107028.6 


0.29(2) 


14.3(5) 


J= 


=79/2-77/2 


109808.4(10) 


109809.0 


0.09(3) 


14.5" 


B1394 


N= 


=29-28 a 


80878.5(4) 


80878.6 


0.34(6) 


14.2(5) 


N= 


=29-28 b 


80889.5(4) 


80889.9 


0.34(6) 


15.2(6) 


N= 


=30-29 a 


83667.0(3) 


83667.6 


0.24(3) 


14.5" 


N= 


=30-29 b 


83678.6(3) 


83678.6 


0.34(3) 


14.5" 


N= 


=31-30 a 


86456.6(3) 


86456.6 


0.34(2) 


14.4(3) 


N= 


=31-30 b 


86467.5(3) 


86467.4 


0.42(2) 


14.4(3) 


N= 


=32-31 a 


89245.9(3) 


89245.6 


0.30(2) 


14.8(2) 


N= 


=32-31 b 


89255.9(3) 


89256.1 


0.33(2) 


14.5" 


N= 


=33-32 a 


92034.1(10) 


92034.5 


0.35(4) 


14.5" 


N= 


=33-32 b 


92044.1(10) 


92044.7 


0.41(4) 


14.5" 


N= 


=34-33 a 


94823.4(10) 


94823.5 


0.32(2) 


14.5" 


N= 


=34-33 b 


94832.8(10) 


94833.4 


0.33(2) 


14.5" 


N= 


=35-34 a 


97612.8(20) 


97612.4 


0.29(4) 


14.5" 


N= 


=35-34 b 


97623.8(20) 


97622.0 


0.26(4) 


14.5" 


N= 


=36-35 a 


100401.8(5) 


100401.3 


0.25(2) 


14.5" 


N= 


=36-35 b 


100412.0(5) 


100410.5 


0.27(2) 


14.5" 


N= 


=37-36 a 


103189.6(5) 


103190.2 


0.22(2) 


14.5" 


N= 


=37-36 b 


103197.6(5) 


103199.1 


0.30(2) 


14.5" 


N= 


=38-37 a 


105978.6(8) 


105979.0 


0.25(4) 


14.2(5) 


N= 


=38-37 b 


105986.1(8) 


105987.6 


0.27(4) 


14.9(5) 


N= 


=39-38 a 


108767.3(15) 


108767.8 


0.10(3) 


14.5" 


N= 


=40-39 a 


111557.0(10) 


111556.6 


0.14(3) 


14.5" 




=41-40 a 


114344.9(6) 


114345.4 


0.18(3) 


14.5" 


N= 


=41-40 b 


114353.2(15) 


114352.9 


0.17(5) 


14.5" 



Table 1 — Continued 



B1401 



J= 


=59/2-57/2 


82677.5(8) 


82677.7 


0.20(2) 


14.5" 


J= 


=61/2-59/2 


85479.5(5) 


85479.3 


0.21(2) 


14.5" 


J= 


=63/2-61/2 


88280.0(20) 


88280.8 


0.24(4) 


14.5" 


J= 


=65/2-63/2 


91082.0(2) 


91082.2 


0.33(3) 


14.5" 


J= 


=67/2-65/2 


93883.1(10) 


93883.6 


0.11(3) 


14.5" 


J= 


=69/2-67/2 


96684.2(10) 


96684.7 


0.09(2) 


14.5" 


J= 


=71/2-69/2 


99486.0(10) 


99485.8 


0.11(2) 


14.5" 


J= 


=73/2-71/2 


102287.0(6) 


102286.8 


0.19(3) 


14.7(6) 


J= 


=75/2-73/2 


105087.6(6) 


105087.6 


0.10(2) 


14.5" 


J= 


=77/2-75/2 


107889.0(10) 


107888.3 


0.07(2) 


14.5" 



Note. — Number in parentheses are la uncertainties in units 
of the last digits. A superscript "a" indicates that the linewidth 
parameter Vexp has been fixed. 



Table 2. Derived rotational constants 



Series 


B(MHz) 


D(Hz) 


Nlines 


J-range 


B1389 


1388. 860( 2) 


33(1) 


13 


8, 29-40 


B1390 


1389. 878( 7) 


-35(3) 


9 


59/2-79/2 


B1394 


1394.609(10) 


32(4) 


22 


29-41* 


B1401 


1401.559(26) 


139(7) 


7 


59/2-75/2 



Note. — *) B1394 corresponds to a ■^S vibronic state 
with 7=15.2(13) MHz and 7d=-l. 53(35) KHz. Number in 
parentheses are la uncertainties in units of the last digits. 



